Supersolidity from defect-condensation in the extended boson Hubbard model 
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We study the ground state phase diagram of the hard-core extended boson Hubbard model on the square lattice 
with both nearest- (nn) and next-nearest-neighbor (nnn) hopping and repulsion, using Gutzwiller mean field 
theory and quantum Monte Carlo simulations. We observe the formation of supersolid states with checkerboard, 
striped, and quarter-filled crystal structures, when the system is doped away from commensurate fillings. In the 
striped supersolid phase, a strong anisotropy in the superfluid density is obtained from the simulations; however, 
the transverse component remains finite, indicating a true two-dimensional superflow. We find that upon doping, 
the striped supersolid transitions directly into the supersolid with quarter-filled crystal structure, via a first-order 
stripe melting transition. 

PACS numbers: 75.10Jm,05.30Jp,74.20.-z,67.40.-w,67.80.-s 
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I. INTRODUCTION 

In 1956, Penrose and Onsager 1 first posed the question of 
whether one could expect superfluidity in a solid - a super- 
solid state - with coexisting diagonal and off-diagonal long- 
range order. They showed that for a perfect crystal, where 
the wave function of the particles is localized near each lat- 
tice site, superfluidity does not occur at low temperature. 
Later it was proposed 2 ' 3,4 that fluctuating defects in imper- 
fect crystals can condense to form a superfluid, and a su- 
persolid state (with both superflow and periodic modulation 
in the density) emerges. In 2004, Kim and Chan reported 
signatures of superfluidity in solid 4 He in torsional oscilla- 
tor experiments, 5 where a drop in the resonant period, ob- 
served at around T ~ 0.2K, suggested the existence of a 
non-classical rotational inertia in the crystal. 4 Following the 
discovery by Kim and Chan, many experiments and theories 
have attempted to explain this fascinating observation; the sit- 
uation remains, however, controversial ] 6 ! 7 ! 8 ' 9 ' 10 ' 11 

On the other hand, with improvements of quantum Monte 
Carlo (QMC) methods, the origin of supersolid phases can 
be studied exactly in both continuum and lattice models. Ex- 
otic quantum phases, including supersolids, are highly sought- 
after in lattice models, particularly those that may be realized 
by loading ultra-cold bosonic atoms onto optical lattices i 12,13 
The generation of a Bose-Einstein condensate (BEC) in a 
gas of dipolar atoms 14 with longer-range interactions pro- 
vides one promising route to search for the supersolid state. 15 
The extended boson Hubbard model is the obvious micro- 
scopic Hamiltonian to study these systems, and supersolids 
have been found in this model on various lattices i 16 ' 17 ' 18 ' 19,20 
A simplified phenomenological picture for understanding the 
supersolid phase in these models has been the aforementioned 
"defect-condensation" scenario: starting from a perfect lat- 
tice crystal at commensurate filling, supersolidity arises when 
dopants (particles or holes) condense and contribute a su- 
perflow. In the simplest scenario - hard-core bosons doped 
above commensurate filling - this phenomenological picture 



suggests "microscopic phase separation" between the crystal 
and superfluid sublattices. Recent work on triangular lattice 
supersolids 21,22 has called this simple interpretation into ques- 
tion, since there one apparently finds examples where parti- 
cles on the crystal lattice also take part in the superflow. How- 
ever, frustration complicates the interpretation of these results, 
since the underlying order-by-disorder mechanism facilitates 
supersolid formation at half-filling. 

To more closely study the degree to which defect- 
condensation plays a role in the mechanism behind lattice 
supersolids, we study the extended hard-core boson Hub- 
bard model with both nearest-neighbor (nn) and next-nearest- 
neighbor (nnn) hopping and repulsive interactions on the 
square lattice. The Hamiltonian is 

H = —t 2_, ( a l a j + a i a \) + Vi n i n j — I 1 n i 
-t' ^2 ( a l a j + a i a j) + v 2 ^ nin i' 

{{id)) {{id)) 

where a\ and ai are the boson creation and annihilation oper- 
ators, rii = a|a, is the number operator, and (i, j) denotes the 
nearest- and next-nearest neighbors. The t' = limit 

has been studied previously, and is known to harbor several 
crystal solids, including a checkerboard structure that upon 
doping is unstable towards phase separation. 18,23 A striped 
crystal and a stable striped supersolid were also found in the 
t' = limits 

Here, we study the full Hamiltonian of Eq. ([TJ using both 
Gutzwiller mean field theory (Section|II]i and stochastic series 
expansion (SSE) quantum Monte Carlo (QMC) simulations 
(Section ITTTb based on the directed loop algorithmi^l^ 5 . We 
confirm that the model contains a variety of lattice crystals, in- 
cluding a checkerboard and striped phase at half-filling, plus 
a quarter-filled solid. 26 The nnn-hopping t' is found to sta- 
bilize a checkerboard supersolid away from half-filling; this 
and other supersolid phases with differently broken symme- 
try are studied in detail. In general, we find that although 
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supersolid phases are not stabilized at commensurate fillings, 
they are readily formed upon doping. However, we demon- 
strate that, contrary to the simple phenomenological picture 
of defect-condensation, where the crystal and superfluid sub- 
lattices are clearly distinct, in at least one of our supersolid 
phases particles from the crystal sublattice also participate to 
a large degree in the superfiow. 



H. GUTZWILLER MEAN FIELD APPROXIMATION 

We begin by surveying the ground state phase diagram of 
the model in several limits using Gutzwiller mean-field theory. 
The Gutzwiller variational method 27 is a powerful technique 
for studying strongly correlated system. The ground state of 
an interacting system is constructed from the corresponding 
noninteracting ground state, 
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Here, the f ni are site-dependent variational parameters, which 
can be optimized via minimizing the energy 



En = 
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of the variational state. The \rii) form the local Fock basis at 
site i with particles in state \m) . In the hardcore limit, 
we only need to keep the local states |0) and Physical 
quantities are calculated within the variational ground states. 
In particular, we measure the local density (rij), the density 
structure factor at wave vector q 



(3) 



and superfluidity due to a finite value of (a*). The coexistence 
of superfluid order and Bragg peaks in the structure factor sig- 
nifies supersolidity. 

We begin by studying the phase diagram within the 
Gutzwiller approximation near half-filling in the absence of 
the nnn repulsion (i.e. for Vi = 0). From Fig.Q] we find that 
at the mean-field level, various phases are stabilized already 
by this restricted set of parameters. These include a uniform 
superfluid (SF), a checkerboard solid (cS) with ordering wave 
vector q = (tt, tt) (see Fig.QJl for an illustration), and in par- 
ticular a checkerboard supersolid (cSS) - with coexisting diag- 
onal order and superfluidity, away from half-filling (p, = 2V\). 

As expected, increasing the nn hopping t destroys the so- 
lidity of the system. In particular, the supersolid region, found 
when the system is doped away from half-filling, becomes a 
uniform superfluid at large t. This clearly indicates, that large 
t destabilizes the supersolid state. 

To study the effects of a finite nnn repulsion V 2 , we first 
identify two limiting cases. For the t — V% model, a stable 
supersolid state is the striped supersolid (sSS) 16,17,23 ' 28 with 
ordering wavevectors q = (tt, 0) or (0, tt) if obtained. For the 
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FIG. 1: Mean-field phase diagram of t — t' — Vi model at (a) t = 0, 
(b) t — 0.5, and (c) t — 1.0. See text for discussion of the phases, 
(d) possible quantum solid configurations: checkerboard solid (cS), 
striped solid (sS) and quarter-filled solid (qS). All quantities are in 
units of t' . 



t' — V\ model, a stable supersolid state is the checkerboard 
supersolid. To capture the behavior of the system between 
these limiting regimes, we introduce a parameter < x < 1, 
which interpolates between the two regions, by setting t — 
x, t' = 1 — x, and V2 = 5x (see Ref. 12611 ). In the following, 
we thus work in units of t + 1' = 1. Figure [2] shows Gutzwiller 
mean field phase diagrams for different values of x = 0, 0.5, 
and 1. While for finite V2, half-filling is obtained for /1 = 
2Vi + 2V2, we still take /1 — 2V\ as the abscissa, in order to 
ease a direct comparison to the previous case of V2 = 0. 

The coexistence of both nn and nnn repulsions is expected 
to stabilize various solid states—: at half-filling, with both 
V\ and V2 large, a checkerboard (striped) solid is formed for 
Vi > 2V 2 {Vi < 2V 2 yk we find that for p = 1/4, a quarter- 
filled solid (qS) (shown in Fig.[T]l) emerges. In order to distin- 
guish the different solids, we measure the structure factors at 
reciprocal lattice vectors q = (n, tt), (tt, 0) and (0, tt). In ad- 
dition, for the striped structure where a strong anisotropy due 
to the broken rotational symmetry occurs, the magnitude of 
the difference between O x = S(tt, 0) and O y — S(0, tt) is al- 
most equal to the sum O x + O y . Keeping track of this quantity 
allows us to easily distinguish supersolids with an underlying 
striped crystal, from supersolids with an underlying quarter- 
filled crystal, in which case all three structure factors become 
finite, but the difference \O x — O y \ is zero (see Table|T]for a 
summary). 

The case x = 0, shown in Fig.|2^ corresponds to the t' — V\ 
model, which we already discussed (compare to Fig. la). At 
x = 0.3 (Fig. I2J5), the superfluid phase expands, as the intro- 
duction of t and V2 destabilizes the checkerboard solid. For 
x = 0.7 (Fig. |2};), with the model parameters approaching the 
t — V2 limit, the checkerboard structure disappears. Instead, 
striped and quarter-filled structurs emerge, including a striped 
supersolid (sSS), and a quarter-filled supersolid (qSS)i 2 ^ In 
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FIG. 2: Mean-field phase diagram of t — t' — Vi — V% model at (a) 
x = : t = 0, t' = 1.0, and Vi = 0, (b) x = 0.3 : t = 0.3, t' = 0.7 
and V 2 = 1.5, (c) x = 0.7 : t = 0.7, t' = 0.3 and V 2 = 3.5 and (d) 
x = 1.0 : t = 1.0, t' = and V2 = 5. All quantities are in units of 
t + t'. 
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TABLE I: Order parameters associated with the different crystal or- 
ders. 



the limiting case x = 1 (Fig.|2}J), we find two different tran- 
sition paths from the sSS to the superfluid upon doing. When 
Vi < 2, the striped supersolid enters the superfluid directly. 
When Vi > 2, the qS solid and qSS supersolid regions are 
passed as an intermediate regime, separating the striped su- 
persolid sSS and the superfluid. 

Clearly, these mean-field phase diagrams provide evidence 
for not only the existence of various different supersolid 
phases, but also for the possibility of direct quantum phase 
transitions between them (in particular between the qSS and 
sSS phases). In the next section, using these mean-field results 
as guidance, we turn to quantum Monte Carlo simulations in 
order to study in details the various supersolid phases, as well 
as the transitions between them. 



III. QUANTUM MONTE CARLO RESULTS 

We performed extensive quantum Monte Carlo (QMC) 
simulations of the Hamiltonian Eq. (Q~|i using a variation 
of the stochastic series expansion framework with directed 
loops. 24 ' 25 Correlation functions of density operators are eas- 
ily measured within the QMC, and crystal order is signified 
by peaks in the q-dependent structure factor of Eq. ([3]). The 
superfluid density is measured in the standard way in terms of 



FIG. 3: p„, S(tt, 7t) and the number density p vs. the chemical 
potential p for Vi = 3.5 and L = 16 at an inverse temperature 
j3 — 10. All quantities are in units of i! . 



winding number fluctuations, 



Ml 



(4) 



where a labels the x or y direction and (3 is the inverse temper- 
ature. Typically the stiffness is averaged over both directions, 
unless measured in a striped phase which breaks rotational 
symmetry (as discussed below). In the following, we choose 
(3 large enough to ensure simulation of ground-state proper- 
ties, and the system size is N = L x L. 

We begin by examining the phase transition into the cSS 
state, identified in Fig.Q] In the limit where V% vanishes, Fig. [3] 
shows the behavior of the QMC observables at Vi = 3.5i'. 
For t 3> t'it — 100t') (open symbols), there is a discon- 
tinuity near /i — 2Vi = —51/, where a checkerboard solid 
with finite S(ir, ir) melts into a superfluid with finite p s via a 
first-order transition. This discontinuous jump in the particle 
density near p — 2V\ = —5t' is a clear indication that phase 
separation would occur in a canonical system. 23 In contrast, in 
the limit t <C f(t — O.Olt') (solid symbols), the discontinuity 
disappears, and a smooth decrease in the structure factor as 
holes are doped into the system is accompanied by an increas- 
ing superfluidity. The coexistence of both finite S'(7r,7r) and 
p a in contrast to the case t' = 0, indicates that a checkerboard 
supersolid state is stabilized by the nnn hopping. In order to 
confirm that this is indeed true, we perform simulations with 
finite nn hopping t. In Fig. [4] we show QMC results as a func- 
tion of nn hopping t for various V\ and the chemical potential 
fixed at p, — V\ — — 3.5i'. A supersolid phase emerges for 
V\ > 2.5t', and a checkerboard supersolid to superfluid tran- 
sition occurs as t increases. The smooth nature of the data 
across the transition region suggests that the destabilization of 
the cSS state upon increased t occurs via a continuous phase 
transition. 

Next, we consider the effect of the nnn repulsion V2, as al- 
luded to in Fig. [2] We focus on the results from simulations 
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FIG. 4: S(tv, 7r)(solid symbols) and p 3 (open symbols) vs. the nn 
hopping t. p — 2Vi = —3.5, j3 = 10 and L = 32. The vertical 
dash line divides the two phases: the right-hand side is a superfluid 
state and the left-hand side is a checkerboard supersolid state. All 
quantities are in units of t'. 



FIG. 5: (Color online). The structure factors, p a , and p for x = 0.9, 
V\ — 1.0, L — 16 and j3 = 10. In the striped supersolid (sSS), p s 
shows a strong anisotropy. All quantities are in units of t + 1 . 



performed at x = 0.9, corresponding to t = 0.9, t' = 0.1 
and V2 = 4.5. Three different values of the nn repulsion are 
chosen: V x = 1.0, 3.0 and 4.5. For Vi = 1.0, the domi- 
nant t and V2 render the model close to a t — V2 model, and 
a striped structure is expected (Fig. In Fig. |5j the equiv- 
alence of O x + O y and \O x — O y \ indicates the absence of 
the quarter-filled structure, and indeed, at half-filling, a stable 
striped solid (sS) is formed. Furthermore, upon hole-doping 
away from p = 1/2, a striped supersolid emerges. To assess 
the behavior of the superflow in the sSS, we measured the 
superfluid densities perpendicular and parallel to the actual 
stripe direction. For this purpose, pjr and p\ are defined by 
comparing the magnitude of O x and O y calculated after each 
Monte Carlo step: when O x > O y , the ir-direction winding 
number W x (see Eq. ©) is counted as Wn and W y is counted 
as W±, and vice versa-i. Fig.|5]clearly exhibits a pronounced 
anisotropy of p s in the sSS phase. Upon further hole doping, 
we observe a melting of the crystal structure to a uniform su- 
perfluid (SF). This completes the quantum melting of the sS 
crystal upon doping holes - proceeding to a uniform super- 
fluid state via an intermediate sSS state with coexisting super- 
flow and crystal order. 

To stabilize the quarter-filled solid, we study the system 
with a strong nn repulsion. Fig. [6] shows the results for 
Vi = 4.5. A quarter-filled solid (qS) is stabilized at p = 
1/4, whereas at half-filling, a sS is formed. Doping away 
from quarter-filling with holes, a "quarter-filled" supersolid 
(qSS) state is formed, 29 as signified by the coexistence of the 
quarter-filled crystal structure and superfluidity. Upon fur- 
ther hole-doping, the qSS eventually melts into a SF. Doping 
slightly away from quarter-filling with additional bosons, we 
observe a similar qSS state. With further doping, however, O x 
and Oy, as well as p s , begin to exhibit significant anisotropies. 

Near p — 2Vi ~ —2, the anisotropy is most pronounced, 



and S(ir, tt) vanishes, signifying a sSS state. We thus observe 
two seemingly unique supersolid states with different underly- 
ing crystal structures. A detailed study of the transition region 
in Fig. |7]indicates the presence of discontinuities developing 
in the structure factors and superfluid density at the transi- 
tion. This indicates a first-order phase transition between the 
two supersolid phases, as traversed by varying the chemical 
potential. In a simple phenomenological defect-condensation 
picture, this transition may be interpreted as occurring via the 
first-order melting of one of the crystal sublattices that differ- 
entiate the qSS from the sSS. This interpretation is discussed 
more in the next section. 

In Fig. [8] with slightly smaller V\ = 3.0, the qSS state 
is still observed, yet with a reduced extent; no obvious qS 
crystal is observed at quarter-filling on this lattice size. How- 
ever, the superfluid density, although finite, shows a large 
dip near p — 2V\ — 2.9 where the average particle density 
nears p — 1/4. In order to examine this more precisely, we 
performed simulations at a fixed particle density p = 1/4, 
by carefully adjusting the chemical potential, and restrict- 
ing measurements to those Monte Carlo configurations with 
a particle number that precisely matches p = 1/4. The data 
in Fig. [9] strongly suggests that the superfluid density indeed 
scales to zero in the thermodynamic limit, revealing the ab- 
sence of supersolid behavior at p = 1/4. This observation 
is consistent with the picture of supersolidity in this model 
occurring only away from commensurate crystal fillings, and 
arising due to the superflow of doped defects placed intersti- 
tial to the ordered solid structures. 



IV. DISCUSSION 

Using mean-field theory and quantum Monte Carlo simu- 
lations, we studied in detail the formation of three supersolid 
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FIG. 6: (Color online).The structure factors, p s , and p for x = 0.9 
and Vi = 4.5. The size is L = 16 with inverse temperature = 10. 
Quarter-filled supersolids (qSS) and a solid (qS) are observed near 
p = 1/4. A qSS to sSS transition is observed upon doping (see 
Fig.|7J. p s also shows strong anisotropies inside the sSS phase. All 
quantities are in units of t + 1'. 




FIG. 7: (Color online). Detail of the structure factor difference 
(main) and superfluid density anisotropy (inset) for the simulations 
discussed in Fig. [6] where the discontinuities developing with in- 
creased system size indicate a first-order phase transition between 
the "quarter-filled" and striped supersolids. 



phases, which arise in the hard-core extended boson Hubbard 
model of Eq. (1). For large nnn repulsion, a stable checker- 
board supersolid phase can be observed provided a sufficiently 
strong nnn hopping is present, if the system is doped away 
from commensurate (1/2) filling. As observed in previous 
studies, the nn hopping itself is not sufficient to promote su- 
perflow within the doped checkerboard crystal. The other 1/2- 
filled crystal observed in this model is the striped solid that 
breaks rotational symmetry. Again, upon doping away from 
half-filling, a supersolid state emerges from the striped solid. 



FIG. 8: (Color online). The structure factors, p 3 , and p for x — 0.9 
and V\ = 3.0. The size is L — 16 with inverse temperature is 
f3 = 10. All quantities are in units of t + t'. 
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FIG. 9: Finite size behavior of p s at p = 1 /4 (see Fig.[8j. The solid 
curve is a guide to the eye. 



Furthermore, at lower density and large repulsive interactions 
(Vi and V2), the underlying density order changes to a quarter- 
filled crystal structure in order to avoid the large repulsions. 
At particle density of exactly 1/4, traces of the superflow van- 
ish. 

The above observations lend strong support to the idea of 
a mechanism for supersolidity involving the condensation of 
dopants (defects) outside of the lattice crystal. Indeed, in 
no instance can we successfully stabilize both a finite crys- 
tal order parameter and a superfluid density at any commen- 
surate filling. However, the simple phenomenological pic- 
ture of the doped-defect condensation clearly breaks down 
at least for the striped supersolid phase, where although pjr 
and p\ show a strong anisotropy, pjr remains finite even 
close to the half-filled striped crystal. This demonstrates 
that the superfluidity in the striped supersolid is not merely 
a one-dimensional superflow through one-dimensional chan- 
nels. This finding is similar to observations in other mod- 
els on the square lattice ) 17 ' 18 ' 28 and contrasts to the very weak 
anisotropies observed on a triangular lattice striped supersolid 
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at half-filling^. 

The presence of different supersolid phases in this model 
also raises the interesting possibility of observing direct 
supersolid-supersolid phase transitions. In particular, upon 
tuning x, we studied the intermediate region between the 
t 1 — Vi and the t — V2 model. We find that there is no di- 
rect transition between the checkerboard and the striped super- 
solid orders as x is tuned - there is always a superfluid phase 
present when the repulsions become comparable. 30 This is 
similar to the case at half-filling, where the superfluid emerges 
along the line V\ = 2V 2 without a direct transition between 
the checkerboard and striped solid, even when both V\ and V2 
are large 4^ 

In contrast, we find a direct transition between the qSS and 
sSS states in this model upon tuning [i — 2Vi. A detailed finite 
size study reveals that this supersolid-supersolid phase tran- 
sition is a first-order stripe melting transition. Tuning from 
the sSS towards the qSS by decreasing the chemical poten- 
tial, an abrupt increase in the superfluid density component 
perpendicular to the stripe direction takes place, correspond- 
ing to a jump into the qSS crystal structure. This observa- 
tion lends itself to the interpretation that, upon traversing this 
phase boundary, one of the two occupied sublattices, that con- 
tribute to the striped crystal, abruptly melts into a superfluid 
component, while the other remains its rigidity, and provides 
the underlying qSS crystal structure. It would be interesting 
to compare this mechanism to that observed in a supersolid- 
supersolid phase transition on the triangular lattice, 22 where 
significantly stronger first-order behavior is observed. There 
have also been proposed more exotic mechanisms, where su- 



perfluids transition into non-uniform solid phases at commen- 
surate filling, which may be compared to the current work. 31 

In conclusion, we have found several ground state phases 
of the hard-core extended boson Hubbard model with nn and 
nnn hopping and repulsion on the square lattice. Most notable, 
we find that supersolid states readily emerge when doped 
away from commensurability "near" their associated crystal 
phases, provided sufficient kinetic (hopping) freedom is pro- 
vided. The model thus proves an ideal playground for future 
study of concepts related to doping and the formation of super- 
solidity through the mechanism of condensed defects. Further 
studies are necessary to understand the detailed nature of the 
transitions between these different solid, superfluid and super- 
solid phases, as well as their finite temperature properties. 
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